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LARGE DEVIATIONS AND AUBRY-MATHER MEASURES 
SUPPORTED IN NONHYPERBOLIC CLOSED GEODESICS 
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, Abstract. We obtain a large deviation function for the stationary measures of 

jl^^i twisted Brownian motions associated to the Lagrangians L\ [p, v) = ^gp(v, v) — 

Xu!p{v), where g is a C°° Riemannian metric in a compact surface {M,g) with 
nonpositive curvature, w is a closed 1-form such that the Aubry-Mather mea- 
sure of the Lagrangian L{p,v) = ^gp{v,v) — u)p{v) has support in a unique 
l/i 1 closed geodesic 7; and the curvature is negative at every point of M but at 

the points of 7 where it is zero. We also assume that the Aubry set is equal 
to the Mather set. The large deviation function is of polynomial type, the 
^^ , power of the polynomial function depends on the way the curvature goes to 

LJ ■ zero in a neighborhood of 7. This results has interesting counterparts in one- 

dimensional dynamics with indifferent fixed points and convex billiards with 
flat points in the boundary of the billiard. A previous estimate by N. Anan- 
J^ , tharaman of the large deviation function in terms of the Peierl's barrier of the 

Aubry-Mather measure is crucial for our result. 
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Q ■ 1. Introduction 

10 ' Large deviations of families of measures are important in many physical applica- 
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tions where one would like to estimate observables at exceptional conditions using 
"physical" measures in the phase space. A typical setting in dynamical systems 
would be to estimate an invariant measure supported in a singular set with respect 
to the Lebesgue measure (for instance, a measure minimizing some variational prin- 
^^ , ciple) in terms of a family of absolutely continuous measures containing a sequence 

^ ' converging to this singular measure. There is a vast literature on the subject in 

mathematical physics, assuming in most of the cases a hyperbolicity condition for 
the dynamical system and/or dimension one for the configuration space. The sub- 
ject of the present article is to study large deviations in the non-hyperbolic setting 
and higher dimensions with the help of weak KAM theory (see [I] [2] [2] [1] for re- 
lated results) . We will consider a family of surfaces initially described in [11] which 
have negative curvature everywhere up to the points along a closed geodesic where 
the curvature vanishes. Since the hyperbolicity of orbits arises from non-parallel 
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Jacobi fields, and the curvature vanishes along 7, the orbit corresponding to 7 is 
not hyperbolic. 

As a motivation for studying such problem we point out that there classes of 
problems (which do not present full hyperbolicity) where a special orbit plays 
an important role. For transformations with a fixed indifferent point (like the 
Maneville-Pomeau map, see [35], [331 ) this point is associated with the phenom- 
ena of phase transition and polynomial decay of correlation. For special billiards, 
where a cusp point can make a trajectory stay for arbitrary long time close to this 
point, this is also associated with polynomial decay of correlation [5T. The careful 
analysis of the evolution on time of these special trajectories determines singular 
characteristics of these special cases of dynamics. Here we analyze the phenomena 
of large deviation associated to a special closed geodesic. The family of probabili- 
ties indexed by A which is considered here is associated to critical solutions of the 
Evans action (see [17] [21]), and, the limit, when A — > 00, is usually known as the 
semiclassical limit (see [3] [30]). 

One of the important results we get is the link between Peierls barrier and 
Busemann functions, which allow us to make accurate analytic estimates of the 
former from the geometry of the surface closed to the vanishing curvature geodesic 
7. We point out that our work considers a pathological case, and some kind of 
analytical control over the lack of hyperbolicity is essential to get meaningful large 
deviation estimates. This is similar to the case of the investigation of the ergodic 
properties of dynamical systems with fixed indifferent points [ 35] [33] or billiards 
with cusps [27] . 

We denote by d the distance on the manifold induced by the Riemannian metric. 

From these estimations we get the main result of the paper which is the following: 



Theorem 1. Let {M,g) be a compact surface with K <0 such that: 

(1) There is a closed geodesic 7 where K = whose orbit supports the (unique) 
Aubry-Mather measure of L{p, v) — ^gp{v, v) — LUp{v). The Aubry set and 
the Mather set of the Aubry-Mather measure coincide. 

(2) K < in the complement of 7. 

(3) There exists m > such that for every geodesic /? : (— e, — e) — > M per- 
pendicular to 7 at /3(0) = /3 n 7 we have that m is the least integer where 

A* Aw be the stationary measure for the Aw-twisted Brownian motion, A > 0. Then 3 
D > such that V open ball A C M (which does not intersect the closed geodesic 

7), 

-(l/i?)inf d(x,7)2+f < hm hni^iUA)) 

x£A A— >+oo A 



and 



lim hn{tixUA))<-DMd{x,-,f+' 

A— >+oo A xeA 



To a given form w, corresponds by duality in Mather Theory, a homology class 
[h] (see for instance [T^ [13] [IS])- In our case this [h] is the homology of the curve 
7. A brief account of some basic results of Brownian motion and Aubry-Mather 
theory is made in the first three sections of the paper. 
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Several precise results about large deviations for hyperbolic dynamical systems 
arc known [2] [1]. In the case of a probability supported on a proper invariant 
hyperbolic set, the large deviation function is typically a linear function in the 
distance to the support of the measure; its slope depending on the hyperbolicity 
of the system (this in the context of Lagrangian dynamics on the torus see also 
[1]). The thermodynamical formalism is the main tool used to get large deviation 
formulae in the presence of hyperbolicity. 

Let us give an outline of the proof of Theorem 1. N. Anantharaman in [1] [2] and 
[3] considers the family of measures mentioned in Theorem 1: stationary measures 
for the twisted Brownian motions arising from twisting a Riemannian Laplacian 
by multiples of a closed one form w. By the general theory of harmonic analysis, 
such measures are absolutely continuous probabilities on the configuration space 
(the manifold M). These probabilities approximate the Aubry-Mather measure 
associated to the Lagrangian given by the kinetic energy of a Riemannian metric 
plus the closed form w. Under certain assumptions (uniqueness of the Aubry- 
Mather measure), N. Anantharaman also considers a large deviation principle for 
this family, and exhibits a deviation function which is given by the Peierl's barrier. 
This nice geometric result tells roughly that the deviation function at a point p 
depends on how far from being minimizers of the Lagrangian action are closed 
loops based at p. We consider here an special example of this setting where one can 
have a sharp control of the Peierls barrier in a neighborhood of a certain closed, 
non-hyperbolic geodesic (the curve 7 in Theorem 1). In [4], among other things, it 
is analyzed a similar problem for a hyperbolic periodic trajectory of the geodesic 
flow. General references for the Aubry-Mather theory are [29] [26] [12] [18) . 

Our estimate of the Peierls barrier comes from sharp bounds for the Busemann 
function (see [7]) of the Riemmanian metric (M, g) associated to lifts of the geodesic 
7 in the universal covering. These bounds are obtained by comparing the metric 
{M, g) in a neighborhood of 7 with an annulus of revolution (under the assumptions 
of Theorem 1). Notice that the Peierls barrier is defined in terms of the weak KAM 
solutions of the considered Lagrangian (see Sections 4 and subsequent). Therefore, 
one of the main issues of the proof of Theorem 1 is to relate Busemann functions and 
Peierls barrier. There is a natural generalization of Busemann functions to convex, 
superlinear Lagrangians (see section 4.9 of 12 for instance) in the context of weak 
KAM theory. Such generalized Busemann functions are used to exhibit fixed points 
of the Lax-Oleinik operator (backward and forward) with infinite critical value (in 
[l2] [To] these functions are called Busemann weak KAM solutions). We would like 
to point out that the Busemann functions we use are just the Riemannian ones, we 
do not need to apply this general notion in our argument. 

Let us make some comments about the assumption in Theorem 1 concerning the 
equivalence between the Aubry-Mather set and the Mather set. This condition is 
of topological nature, as observed in [35] for surfaces. Particularly important for us 
are the results of section 4 in [28] where it is considered the case where the Mather 
set is a single periodic orbit: if 7 separates M (see case 1.2) then the two sets 
are equal. In fact, the coincidence of the Aubry-Mather set and the Mather set is 
generic in homology (see Theorem 3 in [28]). 
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Let US finish the Introduction with some further remarks and problems. First of 
all, notice that Theorem 1 describes how a family of absolutely continuous measures 
in the configuration space approaches a Dirac measure concentrated in the so- 
called projected Mather set (see Section 3 for the definition). The terminology 
Mather measure may have two meanings: the projected Mather measure (with 
support on the manifold) and the one in the tangent bundle. So Theorem 1 is 
a large deviation for the projected Mather measure. A more recent stream of 
ideas give us some hints of how we can obtain a L. D. P. for the Mather measure 
on the tangent bundle. The so-called entropy penalized method presented in [22] 
shows other ways to obtain approximations of the Mather measure on the tangent 
bundle by absolutely continuous probabilities. Under some conditions {M — T", 
convex superlinear Lagrangians, uniqueness of the Aubry-Mather measure), a large 
deviation principle for such procedure is described in [23 . Another way to get a L. 
D. P. for the Mather measure on the tangent bundle is via semi-classical limit of 
Wigner functions [24]. We believe that these procedures can also be extended to 
our setting. 

Finally, let us observe that the geodesic flows considered here are particular cases 
of expansive geodesic flows in manifolds with non-positive curvature [15j j5] |31| |32j 
[llj [25) . The topological dynamics of expansive geodesic flows in manifolds with- 
out conjugate points is well understood, it is about the same Anosov topological 
dynamics. Moreover, the universal covering is a Gromov hyperbolic space. How- 
ever, the ergodic theory of expansive geodesic flows with non-positive curvature 
is almost the same ergodic theory of rank one manifolds. The ergodicity of rank 
one manifolds is a very hard open problem, as well as for expansive geodesic flows 
with non-positive curvature. The family of examples considered in Theorem 1 are 
perhaps the simplest non-Anosov geodesic flows of rank one, they are ergodic an 
even Bernoulli by Pesin theory. So we could ask if it is possible to give a sharper 
description of the ergodic properties of invariant measures in this case (Liouville 
measure, Gibbs measures). In [21] the Holder class of the horocycle flow for the 
metrics considered in Theorem 1 is presented. We believe that from these estimates, 
and some of the ideas described in the present paper, one can detect the so called 
concentration of measure phenomena, or even calculate the decay of correlation of 
the Liouville measure. This will be the purpose of a future work. 

The second author thanks the Ecole Normale Superieure de Paris where most 
of this work was done during a sabbatic period of the author. Special thanks to 
Professor Viviane Baladi who made possible the visit of the second author to the 
ENS Paris. Special thanks too to the Mathematics Department of the Universite 
de Nice, where the second author developed part of this work when invited by 
Professor Ludovic Rifford. 

We refer the reader to [H] for general properties of large deviations. 

2. Preliminaries about diffusion and Brownian motion 

2.1. Diff'usion and Brownian motion in a Riemannian manifold. Let (M, g) 

be a compact C°° Riemannian manifold, let {M,g) be its universal covering en- 
dowed with the puUback of g by the covering map. Let A be the Laplace operator 
of {M,g), A be the lift of the Laplace operator to {M,g). 
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The operator A gives rise to a stochastic process in (Af , g), the Brownian motion. 
It is hnked to A in the fohowing way: given x e M, let C{M., M) be the space of 
continuous paths, P^: be the Wiener measure, let Xt : 7 — > j{t) be a realization of 
the Brownian motion starting at x. Then 

P*(/(x))=e^*^(/(x))=E,(/(X,)). 

for every C°° function / : M — > R. The operator P* is the Heat semigroup of 
{M,g), i.e., the solution of the heat equation ^ = jAu in {M,g). 

We refer the reader to [TB] or [33] for general results on Brownian motion and 
diffusions. 

2.2. Tvifisting the Laplacian by closed 1-forms. We refer the reader to p^ for 
general results on twisted Laplacians. 

Let w be a C°° closed 1-form in M. The Laplacian twisted by the 1-forni u! is 

AMi^)) = A,(/>)) 

= e-^p'^A{e^p'^)f{x). 

where p G M is a, base point, ui is any lift of a; to M, / is a C°° function in M and / 
is any lift of /. Twisted Laplace operators are used to study asymptotic properties 
of the number of closed orbits in a fixed homology class of geodesic fiows of negative 
curvature (via Selberg's trace formula) taking i?e(w) — 0. The semigroup 

pt ^ iA^t 

gives the solution of the twisted heat equation S| = \A^u for the twisted La- 
grangian L^(p,u) = \gr,{v,v) -uj{v). 

The twisted Laplacian appears in a natural way when we want to consider the 
Schrodinger operator for a Mechanical Lagrangian to which we add a closed form 
(the magnetic term). 

2.3. Stationary probability for the twisted Brownian motion. The operator 
P^ acts on the space of measures by (P^)*: 



/ fdiP^r^,) = f Plfd^i. 



The action preserves positive measures and there exist A(a;) and a measure n^j such 
that 

(P^)V„=e^("V. 
Let h^{x) = Jj^j K^^{y, x)dfi-i^, where K^^ is the Kernel of P^. 



Theorem 2.1. There exists a (unique up to normalization) measure v^ — f^dx 
that is a fixed point of the twisted Brownian motion 

Q'fix) = e-*'^(-)/i.(x)-ipi(;i_./)(a;). 

The measure i^uj will be called the stationary measure of the twisted Brownian 
motion. 
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3. Preliminaries about Aubry-Mather measures 

3.1. Aubry-Mather measures. Consider M a compact C°° Riemannian mani- 
fold. Let L : TM -> R be a C°° convex, superlinear Lagrangian. The action of L 
in an absolutely continuous curve c : / — ?► M is Al{c) — Jj L{c{t),c'{t))dt. 

Let A4 [L) be the set of invariant probability measures of the E-L flow of L. The 
action of L in M.{L) is defined by 

AL{^i) = I Ldfi. 

The homology class (Mather, Mane) p(/i) of the measure fi is given by 



where w is a closed 1-form. (Recall that the homology group Hi [M, R) is the dual 
of the cohomology group H^{M, R)). 

A measure /i G A4{L) is called minimizing in its homology class if 

ALi^i) = mi{ALiiy), piiy) ^ pifj.)}. 

An Aubry-Mather measure /i is defined by 

AL{fi)^mi{AL{j^),iyeM{L)}. 

The union of the supports of all Aubry-Mather measures is called the Mather 
set for L. 

Theorem 3.1. (Mather, Mane): The support of a minimizing measure is a Lips- 
chitz graph over an invariant set of global minimizers of the action. 

3.2. Critical energy values and minimizing measures. Both globally mini- 
mizing measures in homology and Aubry-Mather measures arise as minimum (and 
hence critical) points of the Lagrangian action on holonomic invariant measures. 
Moreover, the support of an ergodic invariant measure has constant energy, and 
the support of globally minimizing measures are minimizing orbits of the Euler- 
Lagrange flow by Theorem 13.11 So it is natural to expect that the energy levels of 
the supports of Aubry-Mather measures have critical properties somehow. 

Definition 3.1. The critical value c{L) of L (see section 2.1 [12] ) is defined by 

c{L) — sup{v4i+fe(/3) < for some closed curve f3}. 

keWL 

An holonomic invariant measure p is called globally minimizing if A^ (/i) = 
-ciL). 

Definition 3.2. The strict critical value cq{L) of L (see page 798 [13] } is given by 



where (i : Hi{M, 



co(L) = min c(L - a) = -/3(0), 

aeH^{M,R) 



/3{h) = min Al{p). 

p{p.)=h 



Notice that cq{L) > c{L). The strict critical level is the relevant one regarding 
Aubry-Mather measures. 
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Theorem 3.2. 8 ; The support of an Aubry-Mather measure is contained in the 
energy level E — co{L). 

There are many equivalent geometric characterizations of the strict critical level, 
it is, for instance, the infimum of the energy levels containing a globally minimizing 
orbit of the Lagrangian action with nontrivial (real) homology class. This is why 
the strict critical value is the critical value of the lift of the Lagrangian action to 
the abelian cover of the manifold (see [l^, for instance). 

Let us give some examples. The critical value of geodesic flows is clearly 
while the strict critical value is nonzero, if and only if, the first homology group 
of the manifold is nontrivial. The critical value of a mechanical Lagrangian is the 
opposite value of the maximum of the potential, and the Euler-Lagrange flow in 
energy levels above this value can be reparametrized to give the geodesic flow of a 
Riemannian metric (Maupertuis' principle). The strict critical value of L{p,v) = 
^5p(w, v) — ujp{v), where w is a closed 1-form, is ^ || cj \\l, where || to \\s is the stable 
norm. In particular, under our hypothesis, the Aubry-Mather measure of L(p, v) — 
^gp{v,v) — ujp{v) is supported in the closed orbit i'y{t),Y{t)), t G [0,Per(7)], and 
the form uj is dual to the homology class of 7. So the stable norm of w is just the 
period of 7. 

4. Large Deviations and Weak KAM theory 

4.1. Large Deviations of stationary measures and Peierl's barrier. In this 
section we state the main tool we use to obtain the deviation function for the 
stationary measures of Brownian motions twisted by multiples of a close 1-form. 
Let us start with some basic analytic definitions in the context of weak KAM theory. 
Our main references are [T5] [12]. Through the section, M will be a compact C°° 
manifold, and L : TM — > R will be a C°° convex, superlinear Lagrangian. 

Definition 4.1. The LfLX-Oleinik operators T^ , Tf^ . Given a continuous func- 
tion f : M ^ {—oo, 00), and t > 0, define the function Tf^ [f) : M — >■ K by 

Tr(/)(x) =inf[/(7(0))+ / L{^{t)n'{t))dt], 
'< Jo 

where 7 : [0,i] — > M is an absolutely continuous curve with 7(i) = x. Let 

r+(/)(a;)-sup[/(7(t))- / L{^{t),yit))dtl 
1 Jo 

where 7 : [0,i] — ?> M is an absolutely continuous curve with 7(0) = x. 

The Lax-Oleinik operators T^^ form a continuous time semigroup family of op- 
erators, as well as the operators T^ . They enjoy very nice properties (see 19}), 
in particular, the family of T^^ — co{L)t has a fixed point u~ which is a viscosity 
solution of the Hamilton- Jacobi equation (see Definition 7.2.3 and Proposition 7.2.7 
in [18] and also [l2]). 

It is also true that u^ is a Lipschitz function (see Theorem 4.4.6 and Corollary 
4.4.13 [18;), and Lebesgue almost everywhere we have that 

H{x,dxU~) = co(L). 

The function u~ is is differentiable along the projection M.q of the Mather set. 
Analogously, the family of operators T^ + co{L)t has a fixed point u"*" which is 
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Lipschitz and a weak (in the above sense) solution of a Hamilton- Jacobi equation 
(the so called conjugated Hamilton- Jacobi equation) [18]. Let S^ be the set of fixed 
points of Tf~ — co{L)t, and let 5"+ be the set of fixed points of T^ — co{L)t. 

Definition 4.2. A pair of functions u~ S S~ , u^ G S*"*" is called a conjugate pair 
if u^{x) = u'^(x) for every x G M.q. 

According to Fathi [T9j [18], for each weak solution u^ of the Hamilton- Jacobi 
equation there exists a unique solution u+ such that u^ , u+ form a conjugate pair. 
We have that the differences m^ — u^ of conjugate pairs are always nonnegative, 
and they vanish at the projected Mather set. One of the most interesting questions 
in weak KAM theory is wether the set of zeroes of the differences of conjugate 
pairs coincides with the projected Mather set. In this case, we can characterize the 
projected Mather set (and hence the set of global minimizers of the action) as the 
set of true critical points of the difference of two C^-smooth sub-solutions of the 
Hamilton- Jacobi equation 120' . This property is very useful for applications, and 
motivates the key idea of the proof of our main theorem. 

So it looks very tempting to try to characterize analytically the projected Mafie 
set in terms of the differences of conjugate pairs. However, we have to be careful 
in this point. 

Definition 4.3. The second Peierl's barrier is 

P{x,x) = inf{u^(a;) —u^{x)} 
where the infimum is taken over all conjugate pairs u+, u^ . 

Definition 4.4. As in 18 , we denote the projected Aubry set by 7r(S(L)), the in- 
tersection of the set of zeroes of all conjugate pairs of the Hamilton- Jacobi equation. 

The projected Aubry set is the canonical projection of certain orbits of the 
Euler-Lagrange fiow (see [18] for instance), whose union is called the Aubry set 
S(L). Moreover, the canonical projection tt : E(L) — > 7r(S(L)) is a Lipschitz 
homeomorphism (see [TS], this is a version of the well known Mather's graph The- 
orem) . 

Lemma 4.1. The projected Aubry set contains Mq. 

However, the inclusion of Mq in 7r(S(i)) might be strict: this is often the case 
when there exist globally minimizing connections between different non-wandering 
components of the Mather set. The projected Aubry set has the analytic charac- 
terization we would like to have for the Mather set (see [18] , [19] , [2] [3] and section 
3.7 in 521). 

Proposition 4.2. The second Peierl's barrier P{x,x) is zero, if and only if, x is 
in the projected Aubry set. 

Finally, we are able now to state the main result of the section, which is one of 
the main tools used in the proof of Theorem 1. 

Theorem 4.3. [5] [3] Let {M,g) be a compact surface with K <Q, let 7, w be as in 
the assumptions of the main theorem. Then, the measures fi\uj satisfy the following 
large deviation type formulae 

lim -InipxijiA)) < - ini P{x,x), 

A— >-fcxD A x£A 
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for every closed set A e M, and 



lim -ln{fixu{B)) > ~ ini P{x,x), 



for every open set B. 



In other words, we know that the measures pL\t^{A) tend to zero as A — > cxi; 
so, the Peierl's barrier gives an estimate of the logarithmic rate of convergence. 
The logarithmic rates became worst as the sets approach the Aubry set, and, if 
the closure of an open set meets the Aubry set, then the above rates are just 
0. This result can be interpreted in the present situation as a concentration of the 
stationary measures around a Dirac measure in the space of continuous paths which 
assigns measure one to the closed geodesic in the Mather set, and zero to any set 
not containing this geodesic We stated Theorem 14 . 3 1 suited to our purposes, as it is 
in [3]: we are assuming that the Aubry set and the Mather set coincide, and that 
there is a unique Mather measure. A more general result in j2_ grants that only the 
first one of the inequalities in Theorem 14.31 holds. 

5. BUSEMANN FUNCTIONS 

Manifolds with nonpositive curvature are special examples of manifolds without 
conjugate points (the exponential map at every point is nonsingular). So every 
geodesic in M is globally minimizing, and the convexity of the metric yields the 
existence of two lagrangian, invariant foliations whose leaves are locally graphs of 
the canonical projection. A well known way to define such foliations in terms of 
closed 1-forms is through the so-called Busemann functions: given 9 = (p, v) € 
TiM the Busemann function b^ : M — > R associated to 9 is defined by 

b'^{x)= lim {d{x,jg{t))-t), 

where d is the metric on M. From now on we will also denote such distance by d. 
The level sets of b^ are the horospheres He (t) where the parameter t means that 
7e(i) € Hg{t) (notice that ^e{t) intersects each level set of b^ perpendicularly at 
only one point). The next lemma summarizes some basic properties of horospheres 
(which can be found in [31] [S^ [S], for instance). 

Lemma 5.1. Let {M,g) be a compact C°° manifold without conjugate points. 

(1) b^ is a C^ function for every 9. If{M,g) has nonpositive curvature b^ is a 
C^ function for every 9. 

(2) The gradient Wb^ has norm equal to one at every point. 

(3) Every horosphere is a C^^^ , embedded submanifold of dimension n — 1 
^(ji+K jjiQQ^jig K -Lipschitz normal vector field), where K is a constant de- 
pending on curvature bounds. If the curvature of (Af , g) is nonpositive each 
horosphere is a C^ submanifold. 

(4) The orbits of the integral flow of — V6 , Vt ■ ^^ — ^ -^; '^'"^ geodesies 
which are everywhere perpendicular to the horospheres Hg. In particular, 
the geodesic je is an orbit of this flow and we have that 

^HHeis)) = Heis + t) 

for every t, s Cz R. 
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A geodesic /3 is asymptotic to a geodesic 7 in AI if there exists a constant C > 
such that d{f3{t),^{t)) < C for every t > 0. Nonpositive curvature imphes that 
every two integral orbits of — V6^ are asymptotic. Item (2) tells us that Busemann 
functions are special, exact solutions of the Hamilton- Jacobi equation 

where H : T*M — >■ R is just the puUback g of the metric g by the covering map, 
H{p,v) = ■^gp{v,v). Moreover, if {M,g) is a compact surface without conjugate 
points and 76) C M is a lift of a closed geodesic, the set of points x e M where 
b^{x)+b^^{x) = Ois just the set of lifts of closed geodesies homotopic to 7r(70) which 
are axes of Tjg. This elementary observation is crucial for the section: when TT{'yg) 
is unique in its homotopy class the functions 6^, b^^ behave like a pair of conjugate 
solutions of the Lax-Oleinik operator. Namely, b^ + b^^ takes its minimum value 
zero just at the points of je, and it is positive everywhere else. 

6. PEIERL'S BARRIER AND BUSEMANN FUNCTIONS 

Definition 6.1. The Peierl's barrier is the function h : M x M — > M given by 
h{x,v) = lim info,(^CT(x,y){AL+co(L){.a)}, 

where Crix^y) is the set of C^ curves a : [0,T] — > M such that a(0) = x, 
a(T) = y. 

The above definition given by R. Mane [26] is based in an analogous definition 
due to J. Mather [29]. We introduce the Peierl's barrier by two reasons. First of 
all, it is naturally connected to the second Peierl's barrier defined in the previous 
section; and secondly, the Peierl's barrier is defined in a more geometric way than 
the second Peierl's barrier; it's actually really close to Busemann functions. The 
purpose of the section is to describe in detail the major issues. Let us begin with 
some basic properties of the Peierl's barrier (see [IS]). 

Lemma 6.1. Let M be a C°° compact manifold, and L be a C^ Lagrangian that 
is strictly convex and superlinear in each tangent space TpM . Let h{x,y) be the 
Peierl's barrier of L. Then 

(1) The function h{x, y) is Lipschitz continuous. 

(2) A point X is in the projected Aubry set if and only if h{x,x) = 0. 

(3) Given x G M, there exists a sequence 7„ : [0,t„] — >■ M of minimizers of 
the action of L such that 

• 7«(0) = Jnitn) = X for every n> 0. 

• lim„_>ootn = 00. 

• h{x,x) = \imn^oo{AL+co(L){l7i))- 

So h(x, x) vanishes at the projected Aubry set, like the differences of conjugate 
pairs of the Hamilton- Jacobi equation. The geometry of the manifold shows up in 
item (2) of the above lemma: if L{p,v) = ^gp{v,v) — uj, the minimizers are just 
geodesies of (M, g), so item (2) tells us that the value of h{x, x) is the limit of the 
values of the action of L -I- cq {L) evaluated on a sequence of minimizing loops based 
on X whose lengths go to infinity. Looking closer at the relationship between first 
and second Peierl's barriers we have (see for instance [H]): 

Lemma 6.2. The following assertions hold: 
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(1) For every pair of conjugate functions u^ , u'^ we have 

u~ {x) — u'^ {x) < h{x, y) 

for every x,y S M . 

(2) h{x, y) — sup„- „+{u^(y) — u'^{x)} where the supremum runs over all pair 
of conjugate pairs. 

6.1. Static classes and uniqueness of the weak KAM solutions. The dy- 
namics of the set of global niinimizers determines the uniqueness of the solutions 
of the Hamilton- Jacobi equation. In order to be more precise about this assertion 
we need some definitions. Let us recall that a semistatic curve /3 : [a, b] — > M is 
an absolutely continuous curve such that 

^L+c[o](/3) = M/3(a),/3W). 
A static curve a : [a,b] — > M is an absolutely continuous curve such that 

^L+c[o](/3) = -/^(/3(6),/3(a)). 
A static curve is always semistatic. The curve /3 : I — > M is static, if and only 
if, it is static restricted to any interval contained in /. Since h{x,y) -f h(jj,x) is 
the infimum of the action of the Lagrangian at the critical level c[0], such curves 
are minimizers of the action. In the case of mechanical Lagrangians, static curves 
are the niaximum points of the potential, while semistatic curves are projections 
of orbits of the Euler-Lagrange flow which tend to the singularities. Observe that 
[T^ the Aubry set is the set of static curves. So in our case, the only static curve 
is the closed geodesic 7 supporting the Mather measure. 

Two points 61, 02 in TM are in the same static class if 

h{Tr{ei,e2)) + h{n{e2,ei)) = o. 

According to our assumptions, there is only one static class, whose elements are 
the points of 7. We need the following result from [T2] : 

Lemma 6.3. Suppose that the set of static classes is unique. Then there exist (up 
to additive constant) a unique pair of conjugate solutions of the Hamilton- Jacobi 
equation. 

This yields 

Corollary 6.4. Suppose that there exists just one static class. Then h{x, x) = 
P{x,x). 

6.2. Peierl's barrier in terms of the differences between Busemann func- 
tions. Now, we are in shape to prove the main result of the section. We are going 
to link the Busemann functions of the metric g in M to the second Peierl's barrier 
P{x, x). For this purpose we shall prove that Busemann functions are naturally re- 
lated with the Peierl's barrier h{x, a;), and then apply the above results. Throughout 
the subsection, 7 = 7^ will be the geodesic in the statement of Theorem 1, namely, 
the support of the Mather measure, 6 G TiM is the initial condition of 79. We 
choose a lift 7^ of -fg in M, and a tubular neighborhood N of 7^ such that the 
covering map 11 : M — > M restricted to iV is a diffeomorphism into a tubular 
neighborhood A^(7) of 7. 
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Proposition 6.5. Let {M,g) satisfy the assumptions of the main theorem. Then, 

(1) h{x,x) = P{x,x) = b^''{x) + b^{x) for every x ^ N and x ^ N such that 
n(x) = x. 

(2) There exist positive constants A, B such that 

A inf (6-«(p) + b\p)) < hinip),Il{p)) < B inf (6-«(p) + b^'ip)). 

We prove the proposition in many steps. 

Lemma 6.6. Let x G N{'-f), let T > be the minimum period of 7, and let 
j3 : [0,r„] — > N^ be a closed geodesic loop parametrized by arc length such that 

(1) m = a: - /3(r„), 

(2) /3 is homotopic to n['y]\. 

Then there exists a function S{n), with lim„^+cx) S{n) = 0, such that 

\AL+clomil3) - ib-%i) + b\S:))\ < Sin), 

where x Cz N is any lift of x in N . 

Proof Let /3 C N he a lift of /3, and let x — /3{0). Let T^^ be the covering 
translation preserving 7^. Then, the assumption implies that 

/3(t„) = t;^(5). 

Since the curvature of (M, g) is negative but at the points of 7 we have that 

(1) The (unique) geodesic [p, T" (p)] joining p to T"{p) is contained in N (be- 
cause of the convexity of the metric) , 

(2) The minimum distance from /3[0, T„] to 7^ must converge to as n — > +00. 
Namely, given e > 0, there exist n > such that for every p G iV, we have 

inf d(n,jn) < e. 
9e[p,T^",(p)] ^^' '"^ - 

Let fix, +[0, +00) C N he the geodesic asymptotic to 7^ with /3x.+ {0) = x. Let 
fix, -{—00,0] he the geodesic asymptotic to 7^ with fix,-(0) = ^"(.t). Let us 
denote by [p, q] the geodesic joining the points p, q in M. We can assume without 
loss of generality that 7^(0) is the closest point in 7^ to x. Since fix,+ i l^x,- a-re 
asymptotic to 7^ there exist n{n) — > if n — > +00, and a number s„ > satisfying 

dil3x,+ isn), t3x,~i-Sn)) < ^(n). 

Let us consider the broken geodesic q;„ given by 

an = f3x,+ [0,Sn] U [f3x,+ {Sn), l3x,-{-Sn)] U f3x-[-Sn,0]. 

By the convexity of the metric, (5„ is contained in the region bounded by /3[0, T„] , 
7g, and the geodesies [x,7g(0)], [r^^(x),7g(nT)]. 

Moreover, the length 1(a) of a is 2s„ + d{fix,+ {sn), f3x.-{—Sn)), and s„ -^ +00 if 
n — > +00, so by the definition of e we have 

\T„ -2sn\ < e. 

Let Co; > be defined by 

Px,+ iax,+) = HgiO) n f3x,+, 
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and let ax,~ > be defined by 
Let us still define Sn,+ > by 

l3x,+ {Sn) e Hg{Sn. + ), 

and s„(b-<?(£)+;,»(i)),_ > by 

/3x-i-Sn) e H_g{-nT + Sn-). 

Notice that |(sn.+ + *«,-) — nT\ < e. The above definitions and the convexity of 
the metric yield 

Sfi = CLx,-\- ~r ='>n,+ ^ ^x.— ~r Sn.— - 

And observe that 

ax,+ = -h%PxAax.,+)) + h%i), 

«x,- = -b-'{iixA-^x,-)) + b-%T:;^{i)). 

Since by definition, 

(1) h^{PxAax,+)) = 0, 

(2) b-%(ix.-{-ax,-)) = +nT, 

(3) h-'\Tl{i))^h-~\i)-nT, 
we have 

ax,+ + ax,- = {b^\5:) + b\5:)). 
Therefore, the length l{a) ~ 2s„ + d(/32;^+(s„),/3j;^_(— s„)) satisfies 

2s„ = ((02;,+ + aj;,_) + s„^+ + s„,_ 

which yields 

|(2s„-nT)-(6-«^(x) + fe^~(x))| 

= |(s„,+ +s„,_) -nr| < e. 

Hence, the length T„ of /? satisfies 

|(T„-nr)-(6-'^(x) + 6^~(£))|<2e. 

The same estimate holds for the curve /3 that is homotopic to n['-f\. To calculate 
the action of L(p, v) + c[0](L) = i(v, w) — ojp{v) + c[0](L) in /3, take without loss of 
generality || 7 ||st= 1 = T, and observe that 



uj — I uj — n. 

So we have 

^L+c[0](L)(/3) = ^;(/3) - n + \l{fi) = ;(/?) - n. 
Thus, 

I^L+c[o](L)(/3)-((^'-'(S) + 6'>))l<2e. 
As n — > +00, we can take e arbitrarily small, and this implies the Lemma. D 
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Lemma 6.7. There exists a tubular neighborhood V C N{'^) of ^ with the following 
property: Let x ^ V , and let 7„ be the sequence defined in Lemma \6.1\ of closed 
geodesic loops based on x such that h{x,x) = V\uin^oo{-^L+co(L){ln)) ■ Then there 
exists m{x) > such that 7„ C N{'-f) for every n > m{x). 

Proof. This is will follow from the fact that the Aubry set and the Mather set are 
equal to the closed geodesic 7. Indeed, suppose that the statement is not true. 
Then we can choose a sequence xi of points converging to 7 with the following 
property: 

Let 7^ be a sequence of closed geodesic loops based on Xi with h{xi,Xi) = 
lim„^oo(^L+co(L)(7ri)) f Lemma 16. II Then there exist a subsequence 7^., n^ — > 00 
if i — > 00, such that 

(1) 7^. is not contained in A'^(7) for every i, 

(2) lim.^oo \Hxi,X,) - iAL+co(.L)ilnJ)\ = 0. 

Since the loops 7^. are minimizers based at Xi where h{xi,Xi) tends to zero with 
? — ^ 00, and their domains tend to M, there exists a subsequence of them converging 
to a global minimizer /?. The geodesic /3 would have a point outside N^j) and, by 
the continuity of h{x,x), /3 would also contain a point p with h(q,q) = 0. By 
Lemma [6.11 the point q is in the Aubry set which coincides with the Mather set: 
the closed geodesic 7. This is clearly a contradiction. D 

Lemma 6.8. There exists a constant C > such that if x ^ N{--f), then h{x,x) > 
Cmi^^^{b-^{p) + b^ip)}. 

Proof. Since the function h(x, x) is zero just at the points of the Aubry set (the 
closed geodesic 7) and is continuous, this implies that outside the tubular neighbor- 
hood A^(7) it must be strictly positive. In the same way, the function b^^{p) + b^{p) 
is strictly positive outside any tubular neighborhood of 7^. The periodicity of the 
Busemann function of 7^, and the convexity of the metric, imply that the function 
d{p) = b^^ip) + b^ {p) is convex in M and attains its minimum value outside N 
at the boundary of N . This minimum value is positive, so /i(a;, x) and d{x) have 
analogous behaviour and the comparison stated in the lemma follows. D 



Proof of Proposition [ 

By Lemma 16.71 if w3 want to estimate h{x, x) at points x & V <Z N{"f) it is 
enough to consider the loops based at x contained in N{'^), where V is the tubular 
neighborhood of lemma l6.7l Then, lemma [6T6l combined with this observation yields 
item (1). Item (2) follows from Lemma [6781 

Remark 

Proposition 16.51 has some interesting consequences regarding the regularity of 
the Peierl's barrier that is closely related to the regularity of the solutions of the 
Hamilton- Jacobi equation. The proposition tells us that the function h{x,x) is, in 
a tubular neighborhood of 7, as regular as the Busemann functions of any lift of the 
geodesic 7 in M. In particular, if the manifold (M, g) has non-positive curvature, 
the regularity of the Busemann functions and Peierl's barrier is C^ [15] . The local 
regularity of Peierl's barrier and subsolutions of the Hamilton- Jacobi equation close 
to the support of the Mather measure is quite exceptional, it holds for instance when 
the support is a finite collection of closed hyperbolic orbits [6] . 
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7. Proof of the main Theorem 

7.1. Surfaces of revolution as models. The main idea to show Theorem 1 is 
to apply our estimates of the Peierl's barrier in terms of Busemann functions in 
order to obtain exphcit bounds for the deviation function for surfaces of revolution 
with nonpositive curvature. Then, to get bounds for the deviation function for 
the surfaces in the statement of Theorem 1 we use comparison theorems to get 
information about the asymptotic behavior of geodesies. 

We choose polynomial functions as test functions to generate annulus of revolu- 
tion: If (a;,y,z) are the cartesian coordinates in M.^ let us consider the curves 

r{z) = {a + z^+'',0,z), 

where fc € N. (The exponent 2 + fc grants the existence of Gaussian curvature) . We 
would like to point out that 

• If the annulus is C°° then k is even. 

• We can also consider k e IR+, but in this case the annulus might not be 

7.2. Estimates for Busemann functions in surfaces of revolution. To esti- 
mate the Busemann functions we use the Clairaut equation 

r(7(f))cos(0(7(i))) = c 

where 7(f) is a geodesic parametrized by arc length, r{p) is the distance from p 
in the annulus to the revolution axis z , and 9{'j{t)) is the angle formed by the 
geodesic 7 and the parallel z — z('-f(t)) (in our case, c = a). 

It follows from a result by P. Eberlein that the Busemann functions in nonpositive 
curvature are C^, so we can write them as functions of z in the following way (in 
fact, by the symmetries of the annulus, it is enough to consider the generating 
curve): 

f^ 6 

M+(z) + u-(z) = u+(0)+u"(0)+ / <Vu+{t) + Vu-{t),— >dt 

Jo 



/o 
= 2 



5t 

2 / sin(6{t))dt 
Jo 



= 2 / ,r--T:V2a + t2 + kdt 

Jo a + t^+'' 

from which follows the estimate h{z, z) = Cz^+2- + o{z^^^). 

7.3. Comparison theory and estimates for Busemann functions close to 
non-hyperbolic geodesies. The purpose of the subsection is to apply the esti- 
mate of the previous subsection to surfaces of nonpositive curvature (M, g) with 
K <0 like in Theorem 1. Namely, 

(1) There is a closed geodesic 7 of period T where K = Q whose orbit supports 
the Mather measure of L{p, v) = ^gp{v, v) — ujp(v). 

(2) K < in the complement of 7. 
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(3) There exists m > such that for every geodesic /3 : (— e, — e) -^ M per- 
pendicular to 7 at /3(0) = /3 n 7 we have that m is the least integer where 

Such surfaces might not have an annulus of revolution containing the geodesic 
7 as a waist, so the estimates in the previous subsection might not apply immedi- 
ately. We shall show that in fact, the sum of Busemann functions u+(z) + u^(z) 
is a function of the angle formed by the two geodesies through z G M which are 
respectively, forward and backward asymptotic to a lift of 7. Then, using CAT 
comparison theorems wc shall be able to compare this function with its counter- 
part in the model surface considered in the previous subsection. This will yield a 
comparison of the sums of Busemann functions in both surfaces. 

To begin with, let us consider a Fermi coordinate system $ : S"^ x (— e, e) -^ M, 
where S^ is a circle of length T parametrized by arc length, such that $(i, 0) = 7(t) 
for every t € [0, T]. The curves ^t{s) = <&(t, s), s e (— e, e), are geodesies which are 
perpendicular to 7. 

Let Ah be the annulus of revolution generated by rotating the curve 

rh{z) = (a + 6z2+'',0,z), 

around the vertical axis in M^. The meridians of the annulus of revolution are 
geodesies, and we can use the same map $ to parametrize a subset of the annulus 
containing the waist 70 (i) ~ Rt{rh{0)), where Rt is the rotation around the z-axis 
of angle t (We can take a > such that the length of 70 is T). Item (3) above 
implies that 

Lemma 7.1. Given a > there exist a tubular neighborhood Vy and constants 
&i < 62 with 1 62 ^ fell < cr, such that 

where Ki, is the Gaussian curvature of Ai,. 

Proof. From the assumptions on the surface {M,g), the curvature K when re- 
stricted to a geodesic /? : (— e, — e) — > M perpendicular to 7 at /?(0) = /3 n 7 
satisfies 

where m is the least integer fc > where the fc-th derivative of K{P{s)) is different 
from 0. Also from the assumptions, the number m is the same for every geodesic 
/? as before. So we can find two functions of the form /i(s) = a + 618^+™, /2(s) ~ 
a + 625^+", such that 

(1) The annulus of revolution At,-^, A^^ generated by the graphs of /i, /2 re- 
spectively, have waists 71, 72 of period 27ra — T, where T is the minimum 
positive period of 7. 

(2) -^/2(0) < i^(/3(0)) < -^/i(0) < 0, for every geodesic /3 : (-e,e) - 
M perpendicular to 7. 

By the theory of surfaces of revolution, the curvatures of A^. at their waists 7^ are 
just the opposites of the second derivatives of the fidXt = 0. By continuity, there 
exists a tubular neighborhood of 7 in the parametrization $ satisfying the assertion 
of the lemma. D 
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The estimates for Buscmann functions of At are completely analogous (up to 
multiplication by constants) to the estimates for Ai showed in the previous sub- 
section. Now, we state a version of the comparison theorem for angles of geodesic 
triangles, suited for our purposes. It is based in the well known comparison theo- 
rems for geodesic triangles due to A. D. Alexandrov and V. A. Toponogov, which 
compare the geometry of geodesic triangles in a Riemannian manifold (Af , g) with 
the geometry of geodesic triangles in spaces of constant curvature (the book of 
Cheeger and Ebin [3] is a great reference for the subject). We shall use the surfaces 
of revolution Ah described in the previous subsection instead of manifolds with con- 
stant curvature as comparison spaces. Since this is not the usual version of CAT 
theorems, we give a proof of the result at the end of the section. 

Theorem 7.2. (Comparison theorem for angles) Let S = R x (— e, e) be a strip, 
and let Si — {S,gi), S2 = (S", 52) be two C°° Riemannian metrics in S such that 

(1) The curvatures Ki, K2 of Si, S2 respectively, satisfy K2{t, s) < Ki{t,s) < 
for every (i, s) G S. 

(2) The curves c{t) = (i,0), i S M, and at (s) = (^7^)7 s G (— e, e), are geodesies 
in Si and S2 for every i G K. 

(3) c{t) has unit speed in Si and S2, crt has unit speed in Si for every t G M, 
and (To is perpendicular to c{t) in Si and 82- 

Let Ai(t,s) be the Si-geodesic triangles whose common vertices are c{0),ao{s),c{t). 
Let [aQ{s) , c{t)]i be the Si-geodesic joining ao{s) andc{t). Letai{t,s) be the Si-angle 
formed by cro(s) and [ao{s),c{t)]i at the point (To(s). Then 

(1) A2{t,s) <z Ai{t,s) for every t,s. 

(2) We have that 02(^,3) < ai{t,s). Moreover, they coincide if and only if 
Ai(i, s) = A2(t, s) and Ki{p) = K2{p) for every p inside Ai{t, s). 

Now, we can apply Theorem 1 7. 2 1 for angles to bound from above and from below 
the distance from 7 to its asymptotes. Namely, let us take /3 like before, with 
/3(0) = 7(io), and we consider /3(a) = $(p), a G [0,e), t ^ io- Let A{p,t) be the 
geodesic triangle in {M,g) with vertices $(p), 7(^0) = ^{Po), l{t) — ^{Pt)', and let 
Ah(p, t) be the geodesic triangle in Ah with the same vertices (let us remind that we 
are parametrizing a tubular neighborhood of 7 and a subannulus of Ab by the map 
$). Then by Theorem^ if Kbi'Pix)) < K{^{x)) for every x in (-e, -e) x [0, T], 
we have that Ai,{p,t) C A{p,t) for every i G M. So if we let t — >■ 00, the same 
property remains true for the ideal triangles in both annuli with vertices ^(to), 

*(p). 

Therefore, the angles formed by /3 and the geodesic 7^ through $(p) that is 
asymptotic to 7 is at least the angle formed by /3 and the geodesic jb,p in Ab 
asymptotic to 7^, whenever Kb{^{x)) < K{^{x)) for every x. 

This application of comparison theory leads us to a comparison between Buse- 
mann functions in the annuli Ab, V^. 

Indeed, let us recall that the formula for the sum of the Busemann functions in 
the previous subsection was 

u+{s) + u-{s) = u+(0) +u-{0)+ < Wu+{p) + Vu"(p), — > dp = 

Jo 5p 

<Vu+{p)+Vu-{p),— > dp. 
&P 
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The same formula essentially holds for the tubular neighborhood V-y in {M,g), 
taking ^ as the Fermi coordinate vector field tangent to the geodesies ^tis) (we 
assume that the pair j'{t), j^ has the canonical orientation). Namely, given t e 
[0,r]wehave 

pS r 

u+($(i, s)) + u-($(t, s)) = J < Vw+($(f,p)) + Vu-($(t, p)), — > dp. 
So we get 

u+{^{t,s))+u-{^{t,s))^ f {cos{e+{t,p)))+cos{e-{t,p))))dt, 



where 9^{t, s) is the angle formed by the forward asymptote of 7 containing $(t, s) 
with the vector ^(<i>(i, s)) (respectively, 6^{t,s) is the angle formed by the back- 
ward asymptote of 7 and ^($(i, s))). 

Let u^ , u^ be the Busemann functions in Ab, and let 9^{t,s), 9^{t,s) be the 
angles formed by the forward (resp. backward) asymptotes of 70 and the vector field 
£ at the point ^{t,s). By Theorem VH] we have that if Kbi^{t,s)) < K{^{t,s)) 
then 

9+it,s)<e+it,s), 

and 

0-{t,s)<0+{t,s). 
Moreover, if K{<^{t, s)) < Ki,i<^{t, s)) then 

e+{t,s)<e+{t,s), 

0+{t,s)<0+{t,s). 

Replacing these inequalities in the above integrals we get lower and upper bounds 
for the sum of the Busemann functions in V^ in terms of the formulae obtained in 
the previous subsection. This finishes the proof of Theorem 1. 

Proof of Theorem [772] 

Roughly speaking, the less curved is the space the more convex is the distance 
between geodesies and narrower are the angles of geodesic triangles. Such features 
of comparison geometry are well known since the famous CAT comparison theo- 
rems. We shall give an elementary proof using Gauss-Bonnet theorem for geodesic 
polygons for the sake of completeness. 

First some general remarks about angles and lengths. Angles are conformal 
invariants, and since two metrics in an open, simply connected set of the plane are 
conformally equivalent we have that the angles with respect to S*!, 5*2 formed by 
two vectors v, w tangent to a point in S are the same. Moreover, the 5'2-length of a 
curve c C S* is at least its S'l-length, because the decrease of the curvature increases 
the norm of vectors and length. 

The Gauss-Bonnet theorem for a geodesic triangle A with inner angles a, 6, c 
tells us that 

{a + b + c)= / K + ir, 

J A 
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and the Gauss-Bonnet theorem for geodesic quadrilaterals D with inner angles 
o, 6, c, d is 



(a + b + c + d) = / K + 2n. 

Ja 



Now, let us consider t = — T, for T > 0, let us denote [(To{s),c{—T)]i ~ 7*, and 
notice that 7* can be parametrized in terms of the parameter t since these geodesies 
are everywhere transversal to the vertical geodesies at ■ So let 7? (t) = jf n at ■ The 
triangles Ai(—T,s) have two sides in common, the geodesies c{t),t € [— T, 0], and 
cro(r), r e [0, s]. The geodesies jf are the other sides of Ai(t, s), and we claim 

Claim 1: 7I is contained in Ai(— T, s). 

For suppose that 7I is not contained in Ai(— T, s). Then there exist to G {—T, 0], 
So G (0, s) such that 

(1) c7t„iso)eYi(^ii, 

(2) the angle formed by ^7! and ^(^(so) is at least the angle formed by ^7^ 
and cri„(so). 

Since "ffi—T) = c{—T) there exists a minimum value ro < to such that 

(1) 7|(t) is not in Ai(— T, s) for every t G (ro,to)7 

(2) 7l(ro)=7i(ro). 

Then, there exists z > s such that the geodesic ^l is tangent to 7I at some point 
7f(ii), where ti G (ro,to)- Let si < s be such that i7tj(si) = 7i(ti)- Let zi > si 
be such that 72(^1) = atAzi). Let us consider the geodesic quadrilaterals Di, \22 
with the following geodesic sides and angles: 

(1) The sides of Di are the Si-geodesics c{t),t G [tij^o], ft^([0,si]), crt„([0,so]) 
and 7i([ti,to])- The inner angles are ai, the angle formed by c(t) and ct^^; 
02, the angle formed by c(t) and (Tt„:, as, the angle formed by at^ and 7f; 
and ai the angle formed by ut^ and 7^. 

(2) The sides of 02 are the S'2-geodesics c{t),t G [ti,io], o'ti([0, 2:1]), crto([0, so]) 
and 7|([ti,io])- The inner angles are 6j, j = 1,2,3,4, where &i is the angle 
formed by c{t) and crtj^; 62, the angle formed by c(i) and CTtj,; 63, the angle 
formed by ut-^ and 7I; and 64 the angle formed by ut^ and 7|. 

Observe that Di C 02, and that Oi = 02 = 61 = 62 = f- So the sum of the inner 
angles of Di is 

TT + (as + 04). 
The sum of the inner angles of 02 is 

7r+(6s + &4). 
Moreover, since 7^ is tangent to 7I at 74 (ii) — crtj(zi), and z > s, we have that 

6s > as- 

By Gauss-Bonnet we get that 

4 4 f f 

V" 5j - V" flj = (6s - as) + {h - a4) = / K2 - Ki. 
1~i ,=1 "'□2 "'□1 
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(64 -ai) = / K2- Ki- (63 - as) 

JU2 JOi 



This implies that 



< K2- Ki. 

JD2 JOi 

Since, by assumption, the tangent vector of 7I at 72(^0) points into Ai(t, s), we 
know that 64 > 04, or 64 — 04 > 0. So we would have that the difference of the 
integrals in the right is positive. However, 

Claim 2: The number J^ K2 — /^ Ki is negative. 

In fact, we have 

K2- I Ki^ f K2- I Ki- f K2. 

'Da ^Di JDi JDi JDs-Di 

So 



K2- Ki< K2~ Ki. 

/Da "'□1 "'Di JDi 

The integral of Ki is calculated in terms of the area form of Si. If we put both 
integrals in terms of dtds, the coordinates of the strip S, we get 



Ki= Ji{t,s)Ki{t,s)dtds, 

Jo, Jd, 



K2^ / J2{t,s)K2{t,s)dtds, 

where Ji{t, s) is the Jacobian of the coordinate change of the S'^-area form. Since 
K2{t,s) < Ki{t,s) the Jacobians satisfy J2{t,s) > Ji{t,s). Thus, /^ K2 — 
/□ Ki < and the claim is proved. 

The Claim contradicts b^ — a^ > 0, which was a consequence of the assuming 
that 7I is not contained in the triangle Ai(— T, s). This shows Claim 1. 

Theorem 17.21 follows from Claim 1 (that is item (1) in the statement), since 
A2(— T, s) C Ai(— T, s) obviously implies that the inner angles of A2(— T, s) are 
bounded above by the inner angles of Ai(— T, s). Gauss-Bonnet implies metric 
rigidity if the two triangles coincide. The argument for T < is completely analo- 
gous. 

References 



[1] N. Anantharaman. On the zero-temperature or vanishing viscosity limit for certain Markov 

processes arising from Lagrangian dynamics, J. Eur. Math. Soc, 6 (2): , (2004), 207-276. 
[2] N. Anantharaman. Counting geodesies which are optimal in homology, Erg. Theo. and Dyn. 

Syst., 23 (2): , (2003), 353-388. 
[3] N. Anantharaman. Entropie et localisation des fonctions propres, ENS-Lyon, (2006), 

|http://www.math.pol3rtechnique.fr/ nalini/HDR.pdf| 
[4] N. Anantharaman, R. Iturriaga, P. Padilla, H. Snchez-Morgado, Physical solutions of the 

Hamilton- J acobi equation, Disc. Contin. Dyn. Syst. Ser. B 5 (3) (2005), 513-528. 
[5] W. Ballman, M. Gromov and V. Schroeder, Manifolds of Non-positive Curvature, , (1985), 

Birkhausser. 



LARGE DEVIATIONS AND NONHYPERBOLIC AUBRY-MATHER MEASURES 21 

[6] P. Bernard, Smooth critical sub-solutions of the Hamilton- J acobi equation, Math. Res. Lett. 

14 (2007), n. 3, 503-511. 
[7] H. Busemann, The geometry of geodesies, (1955), Academic Press, 
[8] M. J. D. Carneiro, On minimizing measures of the action of autonomous Lagrangians, 

Nonlinearity 8 (6), (1995), 1077-1085 
[9] Cheeger, J., Ebin, D. G. Comparison theorems in Riemannian geometry. North-Holland 

Mathematical Library, Vol. 9. North-Holland Publishing Co., Amsterdam-Oxford; American- 

Elsevier Publishing Co., Inc., (1975), New York. 
[10] G. Contreras, Action potential and weak KAM solutions, Calc. Var. Partial Differential 

Equations 13 (4) (2001), 427-458. 
[11] G. Contreras and R. Ruggiero, Non-hyperbolic surfaces having all ideal triangles of finite 

area, Bui. Braz. Math. Soc. 28, (1) (1997), 43-71. 
[12] G Contreras and R. Iturriaga. Global Minimizers of Autonomous Lagrangians, (2004), To 

appear. 
[13] G. Contreras, R. Iturriaga, G.P. Patcrnain and M. Paternain, Lagrangian Graphs, minimiz- 
ing measures and Mane's critical values, Geom. Funct. Anal. Vol. 8 (1998), 788-809. 
[14] A. Dembo and O. Zeitouni, Large Deviations Techniques and Applications, (1998), Springer 

Verlag. 
[15] P. Eberlein, Manifolds of nonpositive curvature. Global differential geometry, 223-258, MAA 

Stud. Math., 27 (1989), Math. Assoc. America, Washington, DC 
[16] R. Durrett. Stochastic Galculus: A practical Introduction, (1996), CRC-Press, 
[17] L. C. Evans, Towards a Quantum Analog of Weak KAM Theory, Commun. in Math. Phys., 

(2004), 244, 311-334. 
[18] A. Fathi. Weak KAM Theorem and Lagrangian Dynamics, (2004), To appear. 
[19] A. Fathi, Solutions KAM faibles conjuguees et barrieres de Peierls, C. R. Acad. Sci. Paris, 

Ser. I 235 (1997), 649-652. 
[20] A. Fathi, A. Siconolfi, Existence of C^ critical subsolutions of the Hamilton- Jacobi equation. 

Invent. Math. 155, n. 2, (2004), 363-388. 
[21] M. Gerber and V. Nitica, Holder exponents of horocycle foliations on surfaces, Ergodic 

Theory Dynam. Systems 19, no. 5 (1999), 1247-1254. 
[22] D. A. Gomes and E. Valdinoci, Entropy Penalization Methods for Hamilton-, J acobi Equa- 
tions, Adv. Math. 215, No. 1 (2007), 94-152. 
[23] D. A. Gomes, A. O. Lopes and J. Mohr, The Mather measure and a Large Deviation 

Principle for the Entropy Penalized Method, To appear in Communications in Contemporary 

Mathematics. 
[24] D. A. Gomes, A. O. Lopes and J. Mohr, Wigner measures and the semi- classical limit to 

the Aubry- Mather measure, (2009), To appear preprint 
[25] A. O. Lopes, V. Rosas and R. Ruggiero, Gohomology and Subcohomology for expansive 

geodesic flows. Discrete and Continous Dynamical Systems, Volume: 17, N. 2 (2007), 403- 

422. 
[26] R. Mane. Generic properties and problems of minimizing measures of Lagrangian systems, 

Nonlinearity, N. 9 (1996), 273-310. 
[27] R. Markarian, Billiards with polynomial decay of correlations, Ergodic Theory Dynam. 

Systems 24 (1) (2004), 177-197 
[28] D. Massart, Aubry sets vs Mather sets in two degrees of freedom, (2008) To appear preprint. 
[29] J. Mather, Action minimizing invariant measures for positive definite Lagrangian systems. 

Math. Z., N 2 (1991), 169-207. 
[30] S. Nonnenmacher, Some open questions in wave chaos, Nonlinearity, 21, no. 8 (2008), T113- 

T121 
[31] R. Ruggiero. Expansive geodesic flows in manifolds with no conjugate points, Ergodic Theory 

Dynam. Systems 17 (1) (1997), 211-225 
[32] R. Ruggiero, Expansive dynamics and hyperbolic geometry, Bui. Braz. Math. Soc. vol. 25, 

n. 2 (1994), 139-172 
[33] O. M. Sarig. Subexponential decay of correlations. Invent. Math. 150 (2002), 629-653. 
[34] D. Stroock, Partial Differential Equations for Probabilists, (2008), Cambridge Press, 
[35] L. S. Young, Recurrence times and rates of mixing, Israel J. Math. 110 (1999), 153-188. 



22 ARTUR O. LOPES AND RAFAEL O. RUGGIERO 

Departamento de Matematica, Universidade Federal do Rio Grande do Sul, Porto 
Alegre, RS, Brazil 

E-mail address: arturoscar.lopesagmail.com 

Departamento de Matematica, Pontifi'cia Universidade Catolica do Rio de Janeiro, 
Rio de Janeiro, RJ, Brazil, 22453-900 
E-mail address: rorrSmat.puc-r io.br 



